ASYMPTOTIC BIFURCATION SOLUTIONS FOR COMPRESSIONS 
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TRANSITION REGION AND BARRELLING TO A CORNER-LIKE 
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Abstract. Buckling and barrelling instabilities in the uniaxial compressions of an elastic rect- 
angle have been studied by many authors under lubricated end conditions. However, in practice it 
is very difficult to realize such conditions due to friction. In the two experiments by Beatty and his 
co-authors, it is found that there is a transition region of the aspect ratio which separates buckling 
and barrelling. Friction, which prevents the lateral movement of the end cross section, might be 
the cause. Here, we study the compressions of a two-dimensional nonlinearly elastic rectangle under 
clamped end conditions. One of the purposes is to show, under this setting in which the lateral 
movement of the end cross section is limited, that there is indeed such a transition region. We 
achieve this by constructing asymptotic solutions of the field equations. By using combined series- 
asymptotic expansions, we derive two decoupled nonlinear ordinary differential equations (ODE's). 
One governs the leading-order axial strain and another governs both the leading-order axial strain 
and shear strain. By phase plane analysis, the axial strain can be obtained from one of the ODE's 
as the axial resultant force varies. Then an eigenvalue problem can be formulated from another 
ODE, which is solved by the WKB method. It is found that when the aspect ratio is relatively large 
there is only a bifurcation to barrelling which leads to a corner-like profile on the lateral boundaries 
of the rectangle. When the aspect ratio is relatively small there are only bifurcation points which 
lead to the buckled profiles. A lower bound of the aspect ratio for barrelling and an upper bound 
for buckling are found, which implies the existence of the above-mentioned transition region. The 
critical buckling loads from our asymptotic solutions are also compared with those obtained from 
the Euler's buckling formula. 

Key words, nonlinear elasticity, bifurcation, corner-like profile, barrelling, buckling 

AMS subject classifications. 74G10, 74G60, 35B32, 34B16 

1. Introduction. Two kinds of instabilities, buckling and barrelling, may occur 
when a two-dimensional nonlinearly elastic rectangle being compressed uniaxially. 
The instability of the uniaxial compression of a rod is an old problem which was 
firstly studied by Euler who gave the Euler's buckling formula to predict the onset of 
buckling of a rod. However BernouUi-Eulcr beam theory does not take into account 
the effect of transverse shear strain, and the Euler's buckling formula is not valid 
for relative thick rods. Recently, a more sophisticated beam theory is derived by 
Russell and White [T] by assuming that the axial displacement is of the second order 
of the transverse displacement. The existence of solutions is established and the post- 
buckling solution is computed numerically. In this paper we consider both buckling 
and barrelling instabilities and present some new analytical post-buckling and post- 
barrelling results for both thin and moderately thick two-dimensional rectangle under 
clamped boundary conditions. 

Theoretical analysis of this type of compression problems has been carried out 
by many authors under lubricated boundary conditions, such as [2l[3l|4l[5] and those 
listed in their references. In Simpson and Spector study barrelling instability 
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of compression of a solid circular cylinder. In [2\ they obtain the condition of critical 
compression ratios at which the cylinder will barrel. In [3], for a specific material a 
detailed relationship between critical compression ratio and modes n are investigated. 
In [4] , both barrelling and buckling instabilities of compression of a two-dimensional 
elastic rectangle are studied by Davies. It is found that the rectangle will buckle or 
barrel under some compression ratios. Buckling and barrelling behavior depends on 
the sign of a parameter which is related to the strain energy function and compression 
ratio. When this parameter is nonnegative, buckling will occur first. A set of sufficient 
conditions that barrelling occurs first is given as a theorem. In [5], Davies studies 
buckling of a square column and barrelling of a circular cylinder. The largest critical 
compression ratio Xbuc that a square column begins to buckle is compared with the 
largest critical compression ratio A^^^ that a cylinder begins to barrel. And then 
some conditions of which instability occurs first are analyzed. 

In the above theoretical analysis, lubricated boundary conditions at two ends are 
used and this condition is also used in many other literatures such as two recent papers 
[6l[7] concerning buckling and barrelling. In [6j, Goriely et. al. study the compression 
of a three-dimensional incompressible cylindrical tube under axial load. They derive 
a new and compact formulation of the bifurcation criterion for both barrelling modes 
and buckling modes. Also, a nonlinear correction is made for the Euler's buckling 
formula. The asymptotic expression of the critical aspect ratio separating instability 
by barrelling from buckling is obtained. Simpson and Spector in [7j prove a rectangular 
elastic rod will buckle under quasistatic compression at two frictionless ends. This 
is very valuable since seldom results can prove "a second bifurcation of solutions 
actually bifurcates from a known solution branch when the known branch becomes 
unstable" (see |7j, p. 1). Another proof of the existence of nontrivial bifurcated branch 
in finite elasticity is given by Healey and Montes-Pizarro [5] and lubricated boundary 
conditions are also used. 

However, in practice and experiments, it is very difficult to realize such lubricated 
boundary conditions, especially when the external force becomes large. Not so many 
experimental works on both buckling and barrelling have been carried out. Two well 
known experiments are done by Beatty and Hook [9] and Beatty and Dadras [TO] . 
In [S], circular rubber bars with different aspect ratios p (ratio of diameter with 
length for cylinders) are compressed at two lubricated ends. It is found that when p < 
pi = 0.216, only buckling can occur. While, when p > p2 = 0.228, only barrelling can 
occur. There is a "transition region in which the nature of the strut behavior becomes 
increasingly ambiguous" (see [S], p. 630). They believe that there should be a limiting 
aspect ratio value in this transition region which separates instability by barrelling 
from buckling and this value is estimated by 0.222 graphically. In another experiment 
of [To], three different elastomers including circular cylindrical bars, rectangular bars 
and thick- walled tubular columns are also compressed at two lubricated ends. In each 
case a similar "transition region" is found and similar results are concluded. In these 
two experiments, their experimental results have good qualitative agreement with the 
theoretical analysis of Beatty [iT] . It should be noted that there exist frictions at two 
flat ends although they are lubricated in experiments. It is also remarked by Beatty 
and Hook in [9] that "In spite of lubrication, the material consistently tended to adhere 
to the end plates of the test machine" ( [9], p. 628) and "We are unable to assess in 
any case the extent to which end effects may have influenced our measurements" ( [9] , 
p. 630). So far there is rare literature from either experimental results or theoretical 
analytical analysis about compressions of bars when friction is taken into account. 
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Also the effect of frictions on instabilities of bars is not known. 

In this paper, due to mathematical difficulties in description of frictions, one of 
our motivations is to analysize the effect of clamped boundary conditions in insta- 
bilities. A key feature shared by a clamped constraint and friction is that they both 
limit the lateral movement of the end cross section. We aim to reveal, under such 
a constraint, that there indeed exists a transition region as observed in experiments. 
As far as we know, it seems that no other work has considered such two or three 
dimensional clamped boundary conditions in stability and instability analysis. The 
clamped boundary condition used by others before is mainly a one-dimensional con- 
dition which only includes zero lateral displacement and zero slope at ends. Here by 
using combined series-asymptotic expansions, we can obtain asymptotically approxi- 
mate boundary conditions for clamped ends. 

Another motivation is related to the Euler's buckling formula which gives critical 
compressive force when a rod begins to buckle under compression. This formula is 
based on the BernouUi-Euler constitutive equation. Here we derive a new model 
equation for the critical stress values where the axial strain is emphasized. Numerical 
results show that we can give an improvement of the first critical stress value. 

The structure of this paper is arranged as follows. In section [2l we present the 
field equations and traction free boundary conditions. Without loss of generality, 
we study the isotropic compressible hyperelastic Murnaghan material. Due to the 
difficulty of nonlinearity, we assume the aspect ratio is small, i.e., we are considering 
a thin or moderately thick rectangle. Then in section [31 after nondimcnsionalization 
of the field equations and traction free boundary conditions, we can find two small 
parameters and one small variable. Thus in section [4] by using the combined series- 
asymptotic expansions two decoupled nonlinear ODE's are obtained. One governs 
the leading order of the axial strain which can be written as a singular ODE system. 
Another governs both the axial strain and shear strain. These two equations can also 
be obtained by the variational principle, as shown in section [5l In section [6l by using 
the combined series-asymptotic expansions, the clamped boundary conditions amount 
to a set of asymptotically approximate boundary conditions, which will be used to 
determine the axial strain and shear strain. In section [71 we make some bifurcation 
analysis. By phase plane analysis, the solution of the axial strain can be obtained 
under the clamped boundary conditions. An asymptotic solution for the shear strain 
can be obtained through the WKB method. Then the condition determining the 
buckling critical stress value is obtained under the clamped boundary conditions. For 
chosen material constants, numerical calculations show that when the aspect ratio 
is relatively large (within our assumption of a small aspect ratio), there will be no 
buckling mode. But there is a bifurcation to the corner-like profile which is a barrelling 
instability. This corner-like profile is related to Willis instability which is described by 
Beatty in fT2l . In FTSfM] , Dai and Wang have made some analysis and pointed out that 
this kind of instability could be caused by the coupling effect of material nonlinearity 
and geometry size of the rectangle. When the aspect ratio is relatively small, only 
buckling modes can be found and the thinner of the rectangle the more of the buckling 
modes. One important difference between our buckling and previous works is that our 
buckling is not a special kind of buckling while buckling in some of the previous works 
such as [Ills] is a Euler-type buckling which is very special. A major finding based 
on the analytical results is that under clamped boundary conditions there is indeed 
a transition region such that when the aspect ratio is larger than a critical value the 
barrelling instability occurs and when it is smaller than another critical value the 
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buckling instability occurs. This reveals that existence of a transition region, which 
is in agreement with the experimental results of Beatty and Hook and Beatty and 
Dadras [TUj- In section [51 we will make our conclusions. 

2. Field equations. We study the deformation of a two-dimensional rectangle 
composed of a compressible hyperelastic material. Let the length of the rectangle be 
I and the thickness be 2a and let {X,Y) and {x,y) denote the Cartesian coordinates 
of a material point in the reference and current configurations respectively. The axial 
and lateral displacements are denoted by 

U{X,Y)^x-X, V{X,Y)^y-Y, (2.1) 

respectively. Then the deformation gradient tensor F is given by 

F = {Ux + l)e^ (g) Ejf + Uye^ (g) Ey + Vxey (g) Ex + (Vy + l)ey (g) Ey, (2.2) 

where Ex , Ey , and e^; , By represent the orthonormal basis in the reference and current 
configurations respectively. Here we have chosen Ex = e^, and Ey = ej^. 

Without loss of generality, we assume this rectangle is composed of a Murnaghan 
material whose strain energy function has the following form 

$ = ^(TrE)2 + Ai(TrE2) + jyi(TrE)(TrE2) + ^(TrE)^ + ^(TrE^), (2.3) 

where E = (F'^F — I)/2 is the Green strain tensor, A and fi are known as the Lame 
constants, vi, and 1^4 are other constitutive constants, Tr is the trace of a tensor. 

The first Piola-Kirchhoff stress tensor S containing terms up to the third-order 
material nonlinearity for an arbitrary strain energy function can be calculated by a 
formula provided in |15j , which is given below: 

Sij — a]iikdkl + 2^jilknmdkldmn + Q^^filknmqpdkldnmdpq + 0(|dij|''), (2.4) 

Ojia-nm and a^juknmqp are elastic moduh defined by 



and I is the identity tensor corresponding to a natural configuration. 

Here we study a static problem, and the field equations (neglecting the body force) 
are given by 

dx + dY ^ '> 
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(A + 2fi)Uxx + I^Uyy + (A + ii)VxY + (A + 2vi + 2u2)UxxVy 
+ (m + i^i + ^^4)(t/yX6fx + VxVyy + 2UxyVx) 

+ (A + 2/i + z/i + ^u4){UxUyy + VxVxx + UyyVy + UyVyy + 2UxyUy) 
+ {X + H + + 2v2 + ]^va){UxVxy + VxyVy) 

+ (3A + 6/^ + 6z^i + 2z^2 + 2v4)UxUxx + Hi = 0, (2.8) 



(A + ^l)UxY + iJ-VxX + (A + 2^)yyy + (A + 21^1 + 2v2)UxVyy 

+ (m + + ]^V4){UxxUy + Vyy + 2UyVxy) 

+ (A + 2/i + :/i + izy4)(C/yf7yy + UxxVx + C/xVxx + VxxVV + 2VxVxy) 
+ (A + /i + 3j^i + 2v2 + ]^Vi){UxUxY + UxyVy) 

+ (3A + 6Ai + 61^1 + 2v2 + 2i^4)VyyYY + iJa = 0, (2.9) 



where Hi and only include third-order nonlinear terms. Here and hereafter the 
lengthy expressions of Hi{i = 1, 2, • • •, 11) arc omitted and their expressions are given 
in the appendix. 

In addition to the above two governing equations, we assume that there is no 
lateral distributive loading on the lateral boundaries. Then the stress components 
YixY and S^y- should vanish on the lateral boundaries. We have the following traction 
free boundary conditions: 



= ^J.UY + fiVx + in + 1^1 + ]^Vi){UxVx + VxVy) 



+ (A + 2^i + vi + ^i^a)(UxUy + UyVy) + ^(A + 2i^i + Vi)UYV^ 



^{2vi + v4)(^-\JlVx + \ulVx + -^/l + -^xV^) 

+ i(A + 2vi + V4)UyVx + (a* + 2i^i + v^UxVxVy 

+ i(A + 2/i + 71^1 + 2v2 + \v4){V\Vy + UyV^) 

+ (4j/i + 2v2 + i^4)t/"xC/y VV + ^(A + 2/i + 2vi + //4)C/^ 

0, at y = ±a, (2.10) 
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^yY = XUx + (A + 2^)Vy + (A + 2i^i + 2v2)UxVy + (^A + v2)Ul 

+ + + ^i^iWyVx + ^(A + 2/i + 1^1 + ]^va){UI + V|) 

+ i(3A + e^i + 61/1 + 2z/2 + 2vi)V^ + (i^i + v2){iUxV^ + C/i) 

i(M + 2^.2 + vaWxUI + C/x^i) + i(A + M + 4i.2)C/iH' 

+ i(A + 2^ + 7i^i + 2v2 + ^i^4){U^Vy + VIVy) + \{2v^ + i^4)t/y VxVV 

+ (^ + 2t^i + i^4)C^A'C^y^x + ^(A + 2/i + 12i^i + 4z^2 + ^v^Vy 
= 0, at r = ±a. (2.11) 

We will study the bilurcations of the nonlinear PDE's ((2^ and ((2?9| under (|2A0)) 
and (|2.1ip and some end conditions. 

For the solution obtained if one of its deformation gradients has a finite jump 
across a curve in the rectangle, this solution should satisfy the jump conditions. These 
conditions require the force balance and the continuity of the deformation. Namely, 

E+n = E"n, F+1 = F"l, (2.12) 

must hold on the curve of discontinuity. and F^ denote its limiting values at 
a point on the curve of discontinuity and n is the unit normal and 1 are all vectors 
tangent to the curve. 

3. Transformed dimensionless equations. It is difficult to study the non- 
linear PDE's (|2.8p and p.9p together with traction free boundary conditions (|2.10p 
and (|2.1ip . But we can use the combined series-asymptotic expansion method to deal 
with this complicated system. This approach adapted here is similar to those used 
in [T4l[T6lfT7l[T8] . First we introduce a new set of dimensionless quantities through the 
following suitable scalings: 



U — hu, V — hv, X = xl, Y — yl, e — ^, v = (3.1) 



where h is the characteristic axial displacement which can be regarded as the reduction 
of the distance between the ends, e will be treated as a small parameter since here 
the deformation is considered to be small. We assume the rectangle is thin such that 
v is small, say, v < 0.07; this implies that a/l < 0.2646, i.e., the aspect ratio is less 
than 0.5292. For simplicity of notation, we will drop the bar from x and y hereafter. 

After the above proper scalings, the original field equations and traction free 
boundary conditions become: 

riQUyy -I- (1 - Tl0)'Vxy + + e{aiVyUyy + aiUyVyy + UlUyyU 

-\-2aiUyUxy + 2a2VxUxy + asVyVxy + asUxVxy + a^VyU XX ~r Ca^UxUxx 

+a2UyVxx + aiVxVxx) + s^H^ = 0, (3.2) 

"^yy + (1 ^ 'r]a)uxy + r^aVxx + e{aiUyUyy + a^VyVyy + aiVyyUx + a2UyyVx 
^-a^iVyUxy + a^UxUxy + 2a2UyVxy + 2aiVxVxy + a2UyUxx + aiVxUxx 
+aiVyVxx + otiUxVxx) + e^Hi 0, (3.3) 
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rioUy + rjoVx + e{aiUyVy + aiUyU^ + a2VyVx + a2UxVx) + £^((ai - ^)uy 

3 3 

+aeUyvl + 2ajUyVyUx + aeUyul + -{ai - l)ulvx + -^{ai ~ + a^VyU^Vx 

3 1 

+ -{ai-l)ulvx+aQUyvl + -{ai~l)vl) ^Q, at y = ±^fi), (3.4) 

Vy + {1- 2T]a)ux + e(^aiM^ + ^a^v^ + a^VyUx + ^a4w| + a2UyVx + ^aiw^) 

+£^{a-julvy + (as - ^)vy + aru^Ux + SawVyUx + agVyul + awul 

+3{a2 ~ rio)uyVyVx + asUyUxVx + aeVyvl + ajUxvl) = 0, at y = (3.5) 

where % = jf2^,Vi = Af27I''?2 = = 2(1+2^ and a,(i = 1,- • -,48) are 

constants related to constitutive constants which are given in the appendix. 

It is also hard to analyze the equations (j3.2p - (l3.5p directly. In the next section 
by using series expansion and asymptotic reduction, we will obtain much simplified 
equations. 

4. Series and asymptotic reduction. In section [31 we have obtained nondi- 
mensionlized field equations and boundary conditions. Since we consider the case 
that 1/ is small and —y/i' < y < y/iy, y is a small variable. We can see that the 
whole problem depends on two small parameters e and v, one small variable y and 
one variable x. Assume that u{x, y) and f (x, y) are sufficiently smooth and then they 
have the following Taylor expansions in the neighborhood of y = 0: 

u{x, y) = uq{x) + y'^U2{x) + y'^Ui{x) + ■ ■ ■ + 5y ■ {ui{x) + y'^u3{x) H ), (4.1) 

v{x,y) = d{vo{x) + y^V2{x) + y'^V4{x) + ■■■)+ y ■ {vi{x) + y'^v^{x) H ), (4.2) 

where (5 is a parameter, which is a measure of deflection of the central axis. Since 
the purpose here is to deduce the critical loads for buckling when there is a buckling 
instability and compare the values with those obtained from the BernouUi-Euler beam 
theory, we consider the case that 6 is small such that 5 — o{y^). 

Substituting (14. 1|) and (14. 2p into the boundary conditions (|3.4|) and p.Sp . we can 
get the following four equations: 

mui + Vo^ox + i^{SrioU3 + rjoV2x) + e{aiuivi + aiuiuox + a2Vivox + a2U02:W0£c 
+i^{3aiU3Vi + AaiU2V2 + 3aiMif;3 + 3aiU3Uoa; + 2aiU2Uix + aiUiU2x + 3a2V3Vox 

+a2U2xV0x + a2U2xV0x + 2q;2W2W1i; + a2Wla;t^l3; + a2VlV2x + a2U0xU2x)) 

3 

+e^(a6UiVi + 2a7UiViUox + aeUiul^ + 2 ("i ~ ^)vfvox + asViUoxVox 

+ ^{ai - l)ul^vox) + 0{u^e^iy, S^) = 0, (4.3) 

2rioU2 + rjavix + v{4.rioU4 + ryofsx) + e{2aiU2Vi + 2aiU2U(ix + "2^1 ^'la; + a.2UoxVix 

+v{AaiU4Vi + 6aiU2W3 + ^oiiUaUqx + 2aiU2U2x + 3a2W3fi2; + a2U2xVix + a2ViV3x 

3 

+a2UQxVy,x)) + e^{2aQU2vl + 4a7U2t'iUoa; + '^Oi&U2ul^ + -(ai - l)tilvix + a^ViUoxVix 
+ ^(ai - l)ul^v^x) + 0{y\e^v, 5^) - 0, (4.4) 
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wi + (1 - 2r]o)uox + i^(3u3 + (1 - 2rjo)u2x) + e(^a5«? + a4ViUQx + ]^aiul^ 

+v{2aiu\ + 3a5WiU3 + 3a4i;3Uoi; + aiViU2x + a.iUoxU2x + 2a2U2Vix + -j^aivlx)) 

3 11 

+e^(aiot;i + -{a^ - 1 + 2rio)vfuox + {^4- 'no)viulx + ^i'^i ^ 1 + 2?7o)woj,) 

2v2 + (1 - 2r/o)Mix + v{Av4 + (1 - 2770)^32;) + £(2aiUiU2 + 2azViV2 

+ 2Q2U2W0a; + a2UiVix + aiWOx'^ls 

+i/(6aiit2U3 + AaiUiUi + 6q!5W2W3 + Aa^viv^ + Aa4V4Uox + 3a4V3Uix 

+ 2a4V2U2x + a4Wl2:W2a; + Q;4WlU3a; + a4U0xU3x + 4q!2U4W0x + 3q;2U31'12; 
+ 2Q;2M2W2a; + ai?^l2;W2£c + a2UiV3x + aiVQxVOx)) 

+s^{4:aQUiU2Vi + 6aiQvlv2 + AarUiU2Uox + 6(q;4 - 1 + 277o)wiW2Woa; 

3 

+ (-1 + 2r]o + 2a4)v2ul^ + -{ai - 1 + 2r]Q)vfuix + (-1 + 2770 + 2a4)viUoxUix 
3 

+ -(«! - 1 + 2T]a)ul^uix + 6(ai - l)u2ViVox + 2asU2UoxVox 

+3{ai - l)uivivix + asUiUoxVix + 2061^1 woajWix + 2aYUoxVox) 

^-0(^/^e2^/,(52) = 0. (4.6) 

In these equations, if O {1/^,8^ 1^,6^) terms are omitted, there are ten unknowns uq, ■ ■ 
•, M4, vo, ■ ■ •, V4. To have a closed system, we need another six equations. 

Substituting (|4.ip and (|4.2p into the field equation p. 21) , the left hand side becomes 
a series in y. All the coefficients of y^{n = 0,1,2,3,- • •) should be zero and as a 
result we have a set of infinitely-many equations. It turns out that the first three 
equations(n = 0, 1,2) contain only the above-mentioned 10 unknowns by neglecting 
proper higher-order terms. These equation are 

2r]oU2 + (1 - 'llo)vix + UQxx + £(2q;iM2Wi -I- 2Q;iU2M0a; + asViVix + asUOx^lx 

+a4ViUoxx + azUoxUQxx) + e^H^ + 0{5'^) = 0, (4.7) 



6r?oW3 + 2(1 - r]Q)v2x + uixx + £(6aiM3t;i -I- 8aiU2W2 + laiuivj, + Qaiu^uox 
+&aiU2Uix + AaiUiU2x + 6a2U3i'oa; + ^a2U2xVox + 2(q;2 -t- az)v2Vix 
+ {2a2 + a3)uixVix + 2azViV2x + 2a^UQxV2x + 2a4V2UQxx + cx^uixUqxx 

~\'CX4ViUixx ~t~ ^b'^Ox'^lxx 

+ 2a2U2V0xx + CtlVlxVOxx + Ct2UiVixx 

+a2VoxVixx)+e''Hfi + 0{5'^) = Q, (4.8) 



I2770U4 + 3(1 - 770)^32; + U2xx + e(12Q!i7t47;i + 18aiU2V3 + l2axU4Uox 
-|-10ai?/2W2£c + 3(2a2 -I- a3)7;37;i2; + (4q;2 + az)u2xVix + 3Q;37;it;3a; 

+^a'iUQxVix + 3a4W3U0x£c + abU2x'U0xx + Q;4l'lM22:a; + Ol<sU0xU2xx 

+2a2U2Vixx + aivixvixx) + e^Hr + OiS"^) = 0. (4.9) 
Similarly, substituting (14. ip and (14. 2p into another field equation p.3p and letting the 
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coefRcients of y'^, and y^ be zero, we have 

2v2 + (1 - rio)uix + VoVqxx + s{2aiUiU2 + 2a5ViV2 + 2aiV2Uo^ + a^ViUi^ 
+a3UQxUix + 2a2U2Vnx + 2a2UiVix + 2aiVQxVix + a2Ui'«o XX ~r 01\Vqx'^Qxx 
+aiViVQxx + aiUoxVQxx) + £^Hs + 0(5^) = 0, (4.10) 

6^3 + 2(1 - rio)u2x + Vavixx + £{^ol\u\ + ^a^vxVz + ^a^vzu^x + 2azVYU2x 

+ 2azUi^xU2x + 6Q;2U2f;ia; + 2axv\^ + 2Q;2W2U0a;a; + OLxV^xUVlxx) 

+e^HQ + 0(^2) = 0, (4.11) 
12i;4 + 3(1 — 7?o)fi3K + 'nQV2xx + e(18Q;iU2U3 + 12aiUiUi + 18Q;5i;2f3 + 12Q;5'i;iW4 

+ 12a4U4U0a; + 3(03 + 2a4)w3Uia; + 2(2q;3 + a4)v2U2x + 'ia.3UixU2x + 3a3UiU3a; 
+ia3,uoxUzx + l2a2UiVox + 12a2 M3l'la; + 10Q!2U2f2a; + 6Q;ifi2;W2£!; + QOL2U\VZx 
+ &aiV(ixV3.x + 3a2M3M0a;a; + aiV2xUoxx + 2a2U2Ulx;r + OLiVi^Ui^^ + a2MiU2a;£c 
+ Q!lW0a;U2a;a; + ?>aiVzVQxx + aiU2xV0xx + 2aiW2Wla;£!; + + Q;iUlW22,a; 

+«ifo.f2.) + e'^fio + 0(5') =0. (4.12) 

Now the field equations (|3.2p - (|3.3p and boundary conditions (|3.4p - (|3.5p are changed 
into a one-dimensional system of differential equations (|4.3p - (|4.12p for ten unknowns 
Mo,---,M4,wo, •••,i'4 if we neglect O^S"^) in (j47l) - (j4l^ and 0{v'^,e'^v, 6"^) in ((43| - ([4Jl) . 
By using perturbation method, we can solve U2 from (|4.7p . Then substituting it into 
(|4.1ip . we can get W3. Further we substitute U2 and U3 into (|4.9p and we can get U4. 
The results are given below. 

U2 = — r Vix - Uoxx + e(ail'yiWlx + ai2U0KWl2; + aizUOxUOxx + "wWlWOia;) 

2770 2?7o 

+0(£'), (4.13) 

11, ,11 
3 oryo o 6?7c 

+ai8l'lWl£!;a; + ai9U0a;l'l£!;a; + a20VlU0xxx + Q;2lU02;U0£!;a;£!;) + 0((r'^), (4.14) 



W3 = (-77 + 7 )W1X2: + (-7 + 7 )U02:X2: + £("15^13, + OIiqVixUOxx + aiT^toa; 

3 oryo o 6?7o 



.1 In /I In 

M4 = (-77: + 77 jl^l£C2;£c + (-777 + 77^ jUaxxxx + £{a22VlxVlxx + a23W0a;2;Wla;2; 

12 12?7o 24 12?7o 

+a-2AVixUQxxx + a25W0a;2;U0a;2;a; + Oi2&ViVixxx + a27W02;tii2;a;2; + a2SVlU0xxxx 

+a29UoxUoxxxx) + 0{e^). (4-15) 

In these equations O(e^) terms are not needed for the final results and their expressions 
are not written out. 

Similarly, we can obtain V2, U3 and V4 from (|4.10p . (|4.8p and (I4.12p respectively. 
The results are given below. 

/ 1 ^Ox 1 / 

V2 = (-7 + -z-juix - -^Vovoxx + e{-aiUiU2 + asoviuix + asiUoxUix 



1 1 

-a2UiUi 

+a32ViVQxx + assUQxVQxx) + s^Hu + O(e^), (4-16) 



-a2U2Vox - a2UiVix - aiVoxVix - ■^a2UiUoxx - -^aiVoxUQxx 
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)voxxx + 0{e), 



(4.17) 




)xxxx 



+ 0{e) 



(4.18) 



In (I4.17P and (|4.18p . 0(e) terms are not written out as they are not needed in the 
final results. Now we substitute (|4.13p . (|4.14p and (|4.15p into the boundary conditions 
(|4.4p and (|4.5p to obtain two equations in terms of uq and vi if we further omit O(e^). 
And, we substitute ^J5^, (jil^ . (jiTZ)) and into the boundary 

conditions (|4.3p and (14. 6p to obtain two equations in terms of uo,ui,vo and fi if we 
further omit 0{siy). Then we can get four equations with four variables uo, ui, vq and 
vi, which are given below. 

(-1 + 2r,o)v,x - uoxx + Hil + - |)"o....) + e( - a4«i.i. 

x'^lxx ~^ ^^b'^Qxx'^lxx 
+'^Z6VixUqxxx + OI^IUQxxUqxxx + ^SS^^l^'lKKa; + asg^Oicl^licaa; + (^AQViUqxxxx 

+aiiU(ixUoxxxx)) = 0, (4.19) 

+ (1 - 277o)wOa; + J^((^ - riii)vixx + ]^UOxxx) + £(^a5«? + a4WlW0a; 

XX ~r Oi42U0xVlxx 

+ai2VlUQxxx + a43'«0a;M0xa;a;)) = 0, (4.20) 

770U1 + r]oVox + l'((-^ + '>ll)uixx + (y " ?7o)'"02;a;2;) + e{aiUiVi + aiWiMoj; 

3 

+a2ViVox + a2U02;t'0a;) + £^(a6Uii'i + 2ajUiViUox + aQUiul^ + -{ai - l)vlvox 
3 

+ Q!8WlU0a;f0a; + -^{^1 - 1)^0x^0:!:) = 0, (4-21) 

2 2 
-mUlx - VO^Oxx + ^((y ~ y )"l£ca;£c + ("y + y )«02;a;2;a;) + " "2^1^12; 

-a2UoxUix - a2UlWla; - aiWQxl'la: " a2UiU0xx - aiVQxUoxx - aiViVoxx 

3 3 

-aiUoxVoxx) + £^(--(q;i - l)vfuix - asViUoxUix - -{ai - l)ul^uix 

-3{ai - l)uiViVix - asUiUoxVlx - 2Q;6^;lW0a;Wla; " ^a^UOxVOxVlx - asUlVlUQxx 

-3(ai - l)uiUoxUoxx - ^ajvivoxuoxx - ^aeUoxVoxUoxx - ae^fvoxx 
-2a7ViUoxVoxx - cteulx^oxx) = 0. (4.22) 

We find that (I4.19P and (|4.20p are two equations with two variables vi and uq, (I4.2ip 
and ()4.22p are two equations with four variables vi, uq, ui and vq. (Note: in these 
two equations (|4.2ip and (I4.22p we keep terms in order to consider the nonlinear 
effect of the axial strain on the possible deflection.) We also notice that both (|4.19p 
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and (|4.22p can be integrated once. Then we obtain the foUowing two equations: 

(1 - 2?7o)wi + UQrc + J^((-2 + —Phxx + (-^ + y)W0a;2;a;j + £y-^0:4Vi + a^ViUox 
+ 2'^5"0a; + f^(("39 - a35)uoxxVlx - OisgUQxVixx + a40VlU0xxx - -^{an ~ a38)vl^ 
-Q;38Uif;ia;a; - -(q;37 - Oi4i)ul^^ - anUOxUOxxx)^ = A, (4.23) 



VoUl + VoVQx + ^ + yjMlija; + ("^ ^)VOxxx) + e(^a2MlUl + a2WlU0a; 

\ . 3 3 
+aiz;iM02; + aiUoa;foa; j + e^l^ ("i ~ + asuiviuox + -^{ai - l)"!^^ 

+aevlvox + 2ajViUaxVox + aeUQ^VQx) = B, (4.24) 

where A and B are two integration constants. 

By deleting vixx terms, from (|4.20|) and (I4.23P and using a perturbation expansion, 
we obtain 

Vi = (-1 + 2r]o)uox + OI44VUOXXX + e(a45Uoa; + ^{^i-&ula:x + OtAlUOxUOxxx)) , 

(4.25) 

and 

1 A 

uqx ~ 7;i^uoxxx + s{Diul^ + v{-D2ul^^ - 2D2UoxUoxxx)) 77; ^ — ,(4.26) 

3 4(l-7?o)?7o 

where Di and D2 are constants related to constitutive constants and their expressions 
are given in the appendix. 

Similarly, from (|4.2ip and (I4.24p . we obtain 

ui = ~vox + 2v{ri(i - l)vQxxx + e(2(l - rio)uoxVox) + {aiaul^vox) , (4.27) 

and 

1 B 

-7;i^voxxx + suoxVox + e^{Di - l)ul^vox ^ — . (4.28) 

3 4(1 - ?;o)f7o 

In order to find the physical meanings of the two integration constants A and B, 
we consider the axial resultant force T and the shear resultant force Q. The axial 
resultant force T is given by 

T= f ^xxdY. (4.29) 

After expressing T^xx in terms of mq and vi and carrying out the integration, we find 
that 

2ae(A + 2^) = ' + ""^ + ^^^^ 2 + —^""'^^ + 2 + Y ) 

xx'^lx '^39'^0a;^la:a; 
+ a40WlU0a:2:a; - 2(0^34 - a38)^'L " a38?^lWl:Ea; " 2('^37 " a4l)'«0a;a; 



a4lU0a;U0a;a;a;) ) ■ (4.30) 
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Comparing g^Ol) with g^Sl), we see that - A, i.e., ^ = j^t4^^ 

where E — is an elastic modulus. Now we let 7 — and hence 7 is the 

nondimcnsionalized averaged axial resultant force. Then from (j4.26p . we have 



£UOx 



^veuQxxx + £^{piul^ + v{ ~ D2ul^^ - '2,D2UqxUoxxx)^ = 7- (4-31) 



If we use the original dimensional variable by letting = £Uqx, we obtain 

W + DiW^ +a^(^-^Wxx-D2{Wl + 2WWxx)) =7- (4.32) 
Next we consider the shear resultant force Q, which is given by 

Q = ( ^yxdY. (4.33) 



After expressing in terms of vo and ui and carrying out the integration, wc find 
that 

- rioui + rioVQx + v[[ — — + — jui^a; + [-7. 7r)voxxx) 



+e[a2UiVi + a2UiUQx + aiViUQ^ + aiUo^wox 



3 3 

+afivlvo^ + 2a7ViUoxVox + aeUox^o^:)- (4.34) 

Comparing gJl with gJl, we see that 2,e^(A+2M) = ^' i'^-' 5^ = 4(i4o)^o " 

is the non 

force. Then from (|4.28p . we have 



Now we let q = and hence q is the nondimensionalized averaged shear resultant 



-^eSvvQxxx + e'^SuoxVox + e^S{Di - 1)u1^vqx = q. (4.35) 

Again, if we use the original dimensional variable by letting W{X) — euqx and G{X) = 
sSvqx, we have 

{W+{Di-l)W^)G-^Gxx ^q- (4.36) 

Now we finally obtain two decoupled nonlinear ordinary differential equations 
(|4.32|) and (I4.36p . In the next section, we will derive the same two equations by 
using the variational principle. 

5. Euler-Lagrange Equations. In this section, we will derive exact two same 
decoupled equations as ()4.32p and (|4.36p by using the variational principle. From ()2.3p 
we can get the expression of the strain energy $ up to the fourth-order nonlinearity 
which implies that the stress components are up to the third-order nonlinearity and 
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that will conform with our previous derivations: 

^ = + t^Wx + V^) + M^C/y + UyVx + \v^) + XUxVy + 

+(^ + z^i + y2){UlVY + UxV}) + + '^){UxU^ + UxV^ + U^Vy + V^Vy) 

+ (^ + + + y + y)(t^i + V^) + {n + v^+ '^){UxUyVx + UyVxVy) 

+(^^1 + u2){u],Vy + UxV^) + {^ + u, + v^)uW + + i + y + y + ^x) 

+ (2i/i +1/2 + ^-){UxU^Vy + UxV],Vy) + (m + 2j/i + Vi)UxUYVxVY 

+(^ + f + ^ + y + ^-^Wlul + c/iyi + + vM) 

+(^ + + UyVxV^). (5.1) 

The averaged total strain energy per unit length over cross-section is given by 

$ = — / ^dY. (5.2) 
2« J -a ^ ' 

By using the results in the previous sections, we can express ^ in terms of Uo,u\, vq 
and vi , which is given by 



1 2 1 - 2?7o 1 2 / (1 - 2r]o)^ 

8r?o - 8vf' + 4r?o - H"'""'^ ^ - 8r;f + 24(-l + r,o)vf 

-l + 27?o 1 2 , 1 - 27?o 

12(-1 + 770)??^ 4(6-6r?o)r?S 24% - 247?^ 



3-4770 1 3 - 2770 

24(-l + 77o)r/2 24r7o- 247?§ 24(-l + 770)77^ ^ 

/ as 3 Q!4 2 2 '^'5 3 

+n24r;o- 24772 + 8770- 8772 ^'1"°^ + 8r;o - 877^ + 24770 - 2477^ 

^v{divxv\^ + 02Uoxvf^ + 6'37;i7;i 

2;^0a;a; 

+9'tv\vixx + OiViUOxVlxx + ^'g^Ox'^^la^a: + OiqvIuqxxx + ^'lll'lMOxMOraa: + 6'l2Uoa:'"Oxxx)) 
-orlo 1 I2 /4 — 37)0 o 3 — 2770 

[3^"^ + 4^"^"°^ + 8^"°- + "(243^^1^ + -24 + 24r7o ^^"^- 
, 3-2770 ril ril 2 , 1 - 27?o 

+ -24 + 24770 + 2(67?o - 6772)"'"''°"" + 24770 - 247/2^°-- + 24 - 247/0 ' 
1 - 2770 \ , / "^^i 2 , 2 , "2 

^ITa K'a '"OxVOxxx) + £(3 + Z ^"l^^Ox + "j -T^UiViVQx 

24 - 24770 '8?7o-87;2 8770-87^2 4770 - 477^ 

, "2 , "1 2 , Oil 2 ^ , 2/' 22 

4770 - 477g 8770 -877^ 8770-8772 ^6770-1677^ 

Q!7 o 2a6 o o 3(q;i — 1) o 
+ 1 T^Wl^^lWOx + 77 77-T^l"0x - -, , ^UiV^Vo^ 

4770 - 477g 16770-1677^ 8(-l + 770)770 

Q!8 3(Q!i - 1) 2 , 2q:6 2 2 

4% -477^ 8(-l + 770)770 '''' 16770 - 1677^ 
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"7 2 , 2 2 M \ /c: 



47,0 " H 16770 - 1677i^ 

where Oi{i = 1, • • •, 12) are constants related to constitutive constants which are given 
in the appendix. 

By further using ()4.25p and (I4.27|) and omitting higher-order terms, we can reduce 
the above equation to 



+5^\^^vlxx + \euoxvlx + ^(-Di - 

+ ^a'G| + ^WG' + ^{D, - 1)W'G') , (5.4) 



where 

8 ^ 19770 + 107?o^ + 6(-2 + 377o)iJ2 + 2(-6 + 7m)Di 
^ 12-24?;o ' 

^ 8 - 19770 + 107,2 - 67/oL'2 + 2(-6 + 7770)^1 

= 24^0 • ^'-'^ 

The total potential energy per unit area is then given by 

* = / ^dX ~ E [ -fWdX ~ E [ qOdX. (5.7) 
Jo Jo Jo 

So the Lagrangian is given by 

L^^~ E-fW - EqG. (5.8) 
Further by the variational principle, we have the following Euler-Lagrange equations: 
dL d dL d^ dL 

dW" dX dWx ^ 'dX^ dWxx ~ ' , ^ 

(5.9) 

dL d dL 

^ dXdG^ ~ ■ 

If we omit 0{e^S^) and 0{e^6^) in (|5.9|) i . we have the following two equations: 



W + DiW^ - —Wxx - D2a^{Wj. + 2WWxx) = 7, 
o 

{W + {D^ - l)W^)G- yGxx = q. 



(5.10) 



We find that ((57TO)) i and ((57T0)) 9 are exactly (|4?32)) and (|06)) respectively. 

These are two decoupled nonlinear ODE's. This system is called the asymptotic 
normal form equations of the original nonlinear PDE's. After imposing proper end 
conditions, we will study bifurcations of this system. We can study the bifurcations of 
the original complicated nonlinear PDE's through the study of the simpler decoupled 
nonlinear ODE's ((ETU)) . 
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6. Clamped boundary conditions at two ends. In this section, we consider 
the simpUfications of the clamped boundary conditions. As illustrated in Figure \^7H 
two thin parts of the rectangle are clamped by two rigid bodies. The bottom is fixed 
and the top is slidably supported. We will consider the rectangle between two dashed 
lines. Without loss of generality, hereafter we take the length of the rectangle to be 1 
and then 2a is the aspect ratio (ratio of width with length) of the rectangle. Then we 
can propose the following asymptotic boundary conditions for the clamped ends. 
(1) First, for the points A, B, C, D, there is no lateral movement. So we have 

V{0, ±a) = 0, V{1, ±a) = 0. (6.1) 

From the series expression, we have 

V{0,±a) = hv{0,±a) 



= h 



(5wo(0)±a(wi(0) +0^3(0)) + OihSa"^ ,ha^) = 0. (6.2) 



By omitting 0{h5a^, ha^), we have 

voiO) = 0, (6.3) 

and 

vi{0)+a^V3{Q)^Q. (6.4) 
From the perturbation expansions of vi and and (|4.26p . we can reduce (|6.4p to 

13 1 

(- - 30:44)7 + (^2 ^ ^ 3a44){euox) + {Di{-- + 30:44) + a45)(e'"0x)^ = 0. 



Then we have 



(6.5) 

W^(0) = Ai, (6.6) 



where 

3 — 4770 — 6044 ^(3 — 4770 — 6044)2 + 47(1 — 6044) (£'1(1 — 6044) — 2045) 



Ai 



2(2045 -£'1(1 -6044)) 
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is a small negative number depending on the external force 7 and constitutive con- 
stants. We can see that Ai = when 7 = 0. 
Similarly, we can obtain 



i;o(l) = and 1^(1) = Ai. 



(6.7) 



(2) Secondly, for the points A and B, there is no axial movement which leads to 
the following conditions 



f7(0,±a) = 0. 
From the series expansion of [/, we have 
U{0,±a) ^ hu{0,±a) 



uo(0) ± Sa{ui{0) + a^U3{0)) + 0{ha^,h5a^) 



By omitting 0{ha^, hSa^), we have 



and 



uo{0) = 0, 



(6.8) 

(6.9) 

(6.10) 
(6.11) 



From the perturbation expansions of ui and U3 and (j4.28p . we can reduce (j6.1ip to 

G(0) = A2, (6.12) 

where 

A2 = qK, (6.13) 



and K 



-5r)o 



l-(-|+3r7o)Ai-(a-577o+a48 + (-|+5r)o)£»i)Aj ' 

(3) Thirdly, for the points C and D, they should have the same axial displacement, 



I.e., 

U{l,a)^U{l,-a), 
which leads to the following condition 

G(l) = A2. 



(6.14) 



(6.15) 



Now we list the asymptotic end conditions for the asymptotic normal form equa- 
tions (|57TU)) : 



Wo = 0, G = A2, Mo = 0, W = Ai at X = 0, 



uo = 0, G = A2, W^Ai at X ^ 1. 



(6.16) 
(6.17) 



7. Bifurcation analysis. In this section we will make some bifurcation analysis 
of the asymptotic normal form equations (|5.10p . A similar equation as (|5.10P i has 
been derived by Dai and Wang in [131 [13] for the compressions of a cylinder and 
some results are obtained. Here we will use these results to study our problem. It is 
found that there are two different kinds of bifurcation phenomena for chosen material 
constants. One is that there is only bifurcation to the corner-like profile when the 
aspect ratio is relatively large, which is a barrelling instability. Another is that, when 
the aspect ratio is relatively small, only bifurcation to the buckled profile occurs 
instead of barrelling to a corner-like profile. 
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7.1. Barrelling instability: Bifurcation to a corner-like profile. First we 
can rewrite (IS.lOP i as 

'Wx=g, 




Fig. 7.1. 7-iy plot 



9 9 9 




(a) (b) (c) 



Fig. 7.2. Phase planes for different ■y values: (a)7p < 7 < 0; {b)-f = -fp; {c)-fc < 7 < 7p 
-4Sr <7<7c 

The vector field of the system (j7.1l) has a denominator term, which is zero at 
W — —qJ)^, and thus it is a singular ODE system. Equation (|7.ip together with 
(|6.16P a and ((HUT]) 3 form a boundary- value problem of a singular system. In this 
case, the detailed analysis about this singular dynamical system is almost the same 
with that in [T3]. From (|4.26p . it can be seen that the arising of the denominator 
£)2-term is due to the coupling of the material nonlinearity (measured by e) and the 
geometrical size (measured by v). Here we again point out that it is due to the 
interaction between the material nonlinearity and geometrical size which causes the 
bifurcation to the formation of a corner- like profile as noted in [13 and jll]. 

We concentrate on the case of 3D2 > Di > 0, for which it is sufficient to illustrate 
the bifurcation to the corner-like profile. The equilibrium points of the system are 
given by 



g = 0, -y = W + DiW^. 



(7.2) 
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In the W — g phase plane, there is also a vertical singular line W = ^ There are 
three different phase planes as 7 varies, which are shown in Figure [721 where 7c will 
be defined later and 7^ is the value calculated from (j7.2p 9 for W = — i-e. 



1 



1 



6D2 



(7.3) 



The one-dimensional stress strain relationship is given in Figure [7Tl Without loss 
of generality we let a — 0.25, Di — 2.4, D2 = 0.88, the Poisson ratio a — 2{\+fi) ~ 
0.495 to get the graphical results. Then 7^ = —0.103306. For a trajectory in a phase 
plane to be a solution of the boundary value problem of (|7.1|) . (|6.16P zi and ([03)3, a 
necessary and sufficient condition is that it contacts the vertical line W — Ai twice 
and its X-interval is exactly equal to 1. Next, we discuss the solution(s) in each phase 
plane separately. 
Case (a) 7p < 7 < 

In this case there is a unique solution denoted by trajectory 1. This trajectory 
crosses 14^-axis at (WqjO) in the corresponding phase plane. The solution expression 
is given in the following by referring to [19] . 

r 1 1 2 /, , r Wo-Ei \ 



X : 



-{E2 



6D2 



){E2-Ei) 



[{Wo-E2)n 



'{E2 



E2 - El 

m-,)ne,m)) 

/„ Wo -El 
V ' E2-E1' 



0< X < -, 
- 2 



- < X <1, 
2 - - 



(7.4) 



where 11 and F are the elliptic integral of the third kind and the first kind respectively 
and 



El 



-3 - 2D1W0 - v/3(3 + 167i:'i - iPiWo - 4:D'(W^) 



(7.5) 



E. 



-3 - 2DiWo + ^3(3 + - ADiWq - AD'fW^) 



4:D, 



(7.6) 



/3 = 



1 [D^ 



\E2- Ei){W -Wo) 
{Wo- Ei){W - E2)'' 



{E2 + ^)iWo-Ei) 
\ {Wo + ^){E2-Eiy 
(7.7) 



1 1 

2 



As = Ai at X = 0, Wo is determined by 
2 



^ ^{Wo + ^^)iE2-Ei) 



Wo - El 1 
{Wo-E2)Ili9o, " \ m) + {E2 + — )F{eo,m) 
E2 - El 6D2 



(7.8) 



'{E2-Ei){Ai-Wo) 



where On = arcsin 1 , , , , . 

^ (Wo - Ei)iAi - E2) 

Case (b) 7 = 7p 

In this case, there is a trajectory tangent to the singular line in the phase plane 
as shown in Figure [721 (b)- There is a unique solution denoted by trajectory 1 in the 
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phase plane. Trajectory 1 contacts 1^-axis at (Wq, 0). The solution expression is 
given by 



X = { 



-1 



El 



{Wo-E2){f^ 



Wo -El 

Wo_-Ei 

1 

6D2 



o<x<-, 



(7.9) 



where Ei and E2 are the same as ()7.5p and (|7.6p respectively and 



= arcsm < 



^ iWo^Ei)iW+^y 
As = Ai at X = 0, H/q is determined by 
1 1 2(^^0+6^;) 



{^-E2){Wo~Ei) 



2 JiWo-E2)i^-Ei 
where — arcsin 



iWo-E2){ 



Wo -El 



6D2 



El)' 



(7.10) 



(7.11) 



(-to:7--Ei)(Ai-VKo) 



(Wo-Bi)(Ai + 5i^) ■ 

Case (c) 7c < 7 < 7p 

There is a unique solution represented by trajectory 1. It contacts M^-axis at (Wq, 
0). This solution can be expressed in the following form 



X = i 



f3 



w 



Ai 



rr—l^ ^dW, Q<X < -, 

(W -Wo){W^ + sW + t) ' - - 2' 



(7.12) 



1 



V^ + mr 1 
a Jwa V (W^-VFo)(M/2+sVK + t) ' 2 - - ' 



where s 



6-f + 3Wo + 2DiWi 



following equation 



2Di 

1 _ 

2 ~ 



and t 



3 + 2DiWo 
201 ■ 



And Wo is determined by the 




W 



1 

6D2 



{W - Wo){W^ + sW + t) 



dW. 



(7.13) 



Case (d) 7 = 7c 

We define 7c to be the critical stress value such that when 7 = 7c trajectory 2 
in Figure [7721 (c), which starts from A, goes to B, jumps from i? to C and finally 
arrives at D, is also a solution. In this case, there is another solution indexed by 1 in 
Figure [7T2I (c). For trajectory 1, it contacts ly-axis at (Wo, 0), and its solution also 
is expressed by (|7.12l) and Wo is determined by ()7.13p . We find that for this solution 
the value of Wxx a.t X — 0.5 is relatively small. Indeed, for 7 — 7c, Wxx — 1.68199 
at a: = 0.5. 

The solution expression corresponding to trajectory 2 is given by 



W ■- 



E2 + (E2 - El) sinh^ - ( 2arcsinhl/ — + /sf- -x) 

2\ V -E^ +Eo \2 J 



E2 + {E2 - El) sinh^ - ( 2arcsinh 



-El + E2 



-E1 + E2 V 



< X < 



< X < 1, 



(7.14) 
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where 



Di - 9D2 - J3( -Dl + 27Dj + 6DiD2{l + 247D2)) 
T2DJT2 ' ^'-''^ 



Di - 9D2 + J3( -Dl + 27DI + &DiD2{l + 247D2)) 

E2 = . (7.16) 

12D1D2 



As = Ai at X = 0, we have 



1 1/ . , /A1-S2 „ . , -6S7-^2 



which determines the value 7c. For the previously-chosen parameters, we find 7c = 
-0.1038086. 

Alternatively, (|7.14p can be rewritten as 



W = E2 + {E2-Ei) sinh" \ ( 2arcsinhi/ ^2 _! + 1 _ x ) . (7.18) 



7t £^2 



2 \ y —£^1 + i?2 



1 



2 



This is a non-smooth solution whose first-order derivative at X = 0.5 does not 
exist. This non-smooth solution is a weak solution of the singular ODE system ac- 
cording to the definition given in j20J . However, for this solution the jump conditions 
(|2.12p cannot be satisfied exactly and we shall not take it as one solution of the field 
equations. We intend to find the smooth solution which at one point has a very large 
second-order derivative value. 
Case (e) -^ij < 7 < 7c 

When I7I is slightly larger than |7c|, there are two solutions, which both can be 
indexed by 1. Both of them can be expressed by (|7.12p . We denote the crossing points 
with the M^-axis of the two solution trajectories by (M^,0) and (M^,0) respectively 
and Wo and Wq can be determined by (|7.13p . We notice that for one of them the 
point (Wq, 0) is very close to the singular line, which represents a sharp crest profile 
with a large second-order derivative Wxx at X = 0.5. 

For 7 = -0.103809 < 7c = -0.1038086, wc find that Wo = -0.182968 and 
Wg = -0.189393 which is very close to the singular line V = = -0.189394, and 

the corresponding total potential energy values for the two solutions are respectively 

^' = -5.48367 X IQ-^E, * = -5.48513 x lO^^i? (7.19) 

which are obtained from (|5.7|) {G — for a barrelling instability). 

Trajectory 1 with Wq has a smaller energy value and is thus an energetically 
preferred solution. For this solution, we find that Wxx = 115.8 at Z ~ 0.5, indeed 
a very large value. The solution curve corresponding to this trajectory is shown in 
Figure [7731 (a). The surface profile of the cylinder in the current configuration is shown 
in Figure E3l(b). 

By considering the other values of a, we find that when the aspect ratio 0.376 < 
2a < 0.5292 there will be a bifurcation to the corner-like profile. Thus, for the material 
constants chosen here, 0.376 is a lower bound of the aspect ratio for the barrelling 
instability. We also point out that for the barrelling solution obtained, we find there is 
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(a) 



(b) 



Fig. 7.3. (a) W-X plot; (b) Profile of the lateral boundary of the rectangle (x and y are current 
coordinates). 



no nontrivial solution for (|5.10p 9 together with boundary conditions (|6.16P i . (|6.16p 9. 
(|6.17[) i and (|6.17[) 9. This imphes that no buckhng will occur after the barrelling. 

Remark: The solutions obtained in this section satisfy Wm < < (cf. Fig. 
7.ip . i.e., W is in the "hardening" branch. 



7.2. Buckling instability: Bifurcation to a buckled profile. In this sub- 
section, we discuss the bifurcation phenomenon of buckling. Since the solutions for 
W are discussed in detail in the above, in the following discussion, we concentrate 
on equation (I5.10|) 9. We will formulate our eigenvalue problem with the clamped 
boundary conditions. For convenience we rewrite equation ()5.10|) 9 as 



{W + {Di - 1)W^)G ~ —Gxx = q- 



(7.20) 



This equation together with (|6.16P i . (|6.16p 9. (|6.17P i and (|6.17p 9 compose an eigen- 
value problem where 7, contained in 14^, is the eigenvalue and G is the corresponding 
eigenfunction. 

Letting f{X) = W{X) + [Di - l)W{Xf and r = we have 

v3 

t'^Gxx - f{X)G = -q. (7.21) 

Since we set the length I = 1 and we are considering a thin rectangle, r is a small 
parameter. 

We first solve the homogeneous equation 



r^Gxx - f{X)G = 



(7.22) 



by the WKB method, which is a very useful tool for solving bifurcation problems in 
elastic solids (see, e.g. Fu [H]). 

For a small t, the leading-order solution of the homogeneous equation is 



G(X, r) = CiGi(X, r) + C^G^iX, r) (7.23) 
= Ci(-/(X))-3cos r dt + C2{-f{X))-i sin ^lEEEldt^ 



where Ci and C2 are two constants. 
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By the method of variation of parameters, an particular solution (to the leading 
-order) of (|7.2ip is found to be 



Gp{X, r) - ^ r {f{s)f{X))-'^ sin ^^^^J^dtds. (7.24) 



X 



Then the general solution of (|7.2ip . to the leading-order, is given by 



G{X,t) ^Ci{~fiX)ri cos [ ^^^^dt + C2i-fiX))-Um [ ^ ^^^^ dt 

Jo T Jo T 

f"" i f(s)f iX))-i Sin r ^^^^dtds. (7.25) 



T 



X 



By the condition ^J^2, we have Ci = A2(-/(l))3 = qK{-f{l))i. And from the 
condition (|6.16[) i . we have 

X 

G{X, T)dX = e5{vo[x) - t;o(0)) = e5vo{x). (7.26) 
To satisfy the conditions (|6.17P i and (|6.17p 9. we have 



C2A1 + = 0, 

T 



C2B1 + —B2 = 0, 

T 



(7.27) 



where 



Ai = (-/(l))-isin I '^^J^dt, 



A2=KT{cosf ^ ^^^\ t-l)+ I (/(s)/(l))-^sin / ^ ^^^^ dtds, 

Jo T Jo Jl T 

Bi = [\-f{X))-i sin r ^EMdtdX, (7.28) 



B2^KT{-f{l))i I (^f(X))-i cos I ^^^^dtdX 



1 



■J I 

Jo JO 



(/(s)/(X))"4sin / y^^^dtdsdX. 



Jo Jx ' 

To obtain non-trivial solutions, we need the determinant of the coefficient matrix 
of (|7.27p to be zero, i.e., 

A1B2 - A2B1 = 0. (7.29) 

This is the algebraic equation for determining the eigenvalue 7. 

In Figure [73 we plot the curves {A1B2 - A2-Bi)-7 for a = 0.045 and a = 0.06. 
The stress values at intersection points of the curve with 7 axis are the critical stress 
values. There are three critical stress values when a — 0.045 and two when a — 0.06. 
It is expected that more critical stress values will appear when the aspect ratio 2a is 
smaller. 
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-0.08 -0.06 •■^0.04 .--0.02 
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2 
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-0.1 -O.OS" -0.04 -0.02 
(b) 



Fig. 7.4. The {A1B2 - ^2Bl)^ curve for (a) a = 0.045; (b) a = 0.06. 

In Figure 17.51 we give the curve for a ~ 0.25. It can be seen that there is no 
bifurcation point to the buckhng profile before and after the barreUing to the corner- 
hke profile occurs. Actually, when a > 0.0955, there will be no bifurcation to buckling. 
Thus, for the material constants chosen here, 0.191 is the upper bound of the aspect 
ratio for the buckling instability. 



A1B2 -A2B1 
3 



-0.1 -0.08 -0.06 -0.04 -0.02 



Fig. 7.5. The {A1B2 - curve for a = 0.25. 

In Figure FTBl we give the eigen-shapes of the rectangle at bifurcation points when 
a = 0.045. 

Note that the elastic instability of the compression of a bar at two ends was 
first discovered by Euler in 1744 and then by Lagrange and Bernoulli. A historical 
statement about Euler's derivation can be found in [52] by Komkov. According to the 
BcrnouUi-Euler constitutive equation which says that the bending moment is linear 
in the change of curvature (see Antman ;23j), a basic fourth-order linear ODE for the 
transverse displacement of a slender rectangle is obtained, 



Elu^^'^{s) + Au"{s) = 0, 0<s<l, 



(7.30) 



where u is the transverse displacement of the same material point on the rectangle 
between the undeformed configuration and the current configuration, / = J"^ Y^dY = 
is the moment of inertia and A is the absolute value of the compressive applied 
force. 

There is a detailed derivation of (|7.30p in [24j by Landau et. al in the studies of 
small deflections of rod under axial compressions. Here we also can obtain this equa- 
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0.09 




Fig. 7.6. Eigen-shapes of the rectangle obtained by our model equation for a = 0.045; (a) 
Rectangle in the reference configuration; (b) (c) (d) Buckled rectangle at the frist, second and third 
bifurcation points in the current configuration respectively. 



tion without using the BcrnouUi-Euler constitutive equation. First we take derivative 
to (I4.2ip and add this equation to (|4.22p . then we can obtain 

+UlxUox - ViVOxx - Uq^Voxx + UiVix + UiUQxx " VixVQx - UoxxVOx) 

-2uqxVqxUoxx - Vj^vqxx - Uq^vqxx) + {-Vo + 2aio)(wi 

-UOxVOxVlx + UiViUoxx - ViVOxUQxx - VlUOxVOxx)) = 0. (7-31) 

Further the substitution of vi and ui from (I4.25P and ()4.27p leads to the foUowing 
equation if we write in the dimensional variables 

- IVoxxxx + WVoxx + WxVox + (^i - l)(,2WWxVox + W^Voxx) = 0, 

(7.32) 

where W = euox as denoted before and Vq — sSvq which is the transverse displacement 
of the axis of the rectangle. If we assume that is a constant and omit higher-order 
nonlinear terms, we can obtain, 

EIVoxxxx - EWVoxx = 0. (7.33) 

Since A = —EW in the linear case, this equation can be written as 

EIVoxxxx + AVoxjf = 0, (7.34) 

which is the same as (|7.30|) . 

The general solution of (|7.30p is given by 

u{s) = A + Bs + CsmXs + D cos Xs, < s < 1, (7.35) 

A ' 
EI 

boundary conditions 

u(0) = u(l)=0, u'(0) = u'(l) 0, (7.36) 



where A = and A,B,C and D are constants. If we also impose the clamped 

V EI 
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the eigenvalue A should satisfy 

AsinA + 2cosA - 2 0. (7.37) 

Note that the clamped boundary conditions (|7.36p are used by many literatures such 
as Landau and Lifshitz ( [21], p. 72) in which the clamped boundary condition is 
explained as: "The end of the rod is said to be clamped if it cannot move either 
longitudinally or transversely, and moreover its direction (i.e. the direction of the 
tangent to the rod) cannot change". In section [6l we give the asymptotic clamped 
boundary conditions (|6.16p and (|6.17p for a two-dimensional rectangle. There are two 
significant differences between these two clamped boundary conditions. One is that 
our clamped boundary conditions are related to the axial strain W, see (|6.16l) zi and 
(UTTIIs- Another is that the tangent values of the axis at two ends are not zeros and 
is related to constitutive constants and the external forces 7 and g, see (j6.16|) 9 and 
(16.17^ 9. From the derivations of the asymptotic clamped boundary conditions, we can 
see that these differences are caused by the geometry of the rectangle, the effects of 
the nonuniform axial strain and nonlinearity of the problem. 

Under clamped boundary conditions (|7.36p . the nontrivial solutions of (|7.30p are 
given by 

u(s) ^ cf cosAs- 1 + f (sinAs- As)V (7.38) 

V A — sm A / 

where A is the root of (|7.37p and C is constant. 

The smallest positive eigenvalue of (|7.37p is Ai = 2tt. Corresponding to this 
eigenvalue, the external force A is given by 

A = AEIn'^ (7.39) 

which is the Euler's buckling formula with two ends clamped. By this Euler's buckling 
formula, the first critical stress value (i.e. the smallest absolute ciritcal stress value) 

1 A 4 2 2 

IS given by -— ^ = ^tt^o . 

In Figure 17771 we give the eigen-shapes of the rectangle at critical stress values for 
this 4th-order hnear ODE (fTSO)) . 




2 




Fig. 7.7. Eigen-shapes of the buckled rectangle for the first three bifurcation points obtained 
from the Euler 's buckling formula. 

If we compare Figure FTBl with Figure FtTtI wc find that the eigen-shapes of rectangle 
in the buckled state at the first three bifurcation points coincide. Then it is worthy 
to compare their corresponding critical stress values. 

Before this comparison, we first introduce Davies's work [4J. Davis in [3] studies the 
buckling and barrelling instabilities in the compression of a two-dimensional elastic 
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rectangle. In [1], the traction free boundary conditions at the lateral boundaries and 
the lubricated boundary conditions at two ends are used. At certain compression 
ratio, the bar is always unstable. There are two kinds of instabilities: barrelling and 
buckling. Also, numerical computation can give the critical stress values for different 
barrelling or buckling modes. In the following tables, we compare the critical stress 
values at the bifurcation points for three different models with the same aspect ratio 
of the rectangle and the same kind of material used here. 

In Table 17.11 and 17. 2( we list the absolute values of the first three critical stress 
values for three different models. For the Euler's buckling formula and our model 
equations, they have the same corresponding eigen-shapes shown in Figure 17.71 Ac- 
cording to the formula provided by Davies [4], the bifurcation points in (C) provided in 
Table [73] and [721 correspond to buckling instead of barrelling. From the these tables, 
we find that the absolute critical stress values obtained from [3] are always smaller 
than those obtained by the other models. It is expected that the rectangle under 
the lubricated boundary conditions begins to buckle earlier than under the clamped 
boundary conditions. From the above two tables, we also find that the critical stress 
values obtained by our model equation are very close to those obtained from (|7.30p . 
especially for the first critical stress value. 

Table 7.1 

The averaged nondimensionalized force at the first three bifurcation points when a = 0.045 
calculated from (A) Euler's buckling formula, (B) our model eguation and (C) Davies's method. 





1 


2 


3 


(A) |7i| 


0.0266479 


0.054515 


0.106592 


(B) I72I 


0.0268605 


0.0516636 


0.0922735 


(C) |73| 


0.00512409 


0.0170796 


0.0327633 



Table 7.2 

The averaged nondimensionalized force at the first three bifurcation points when a = 0.06 
calculated from (A) Euler's buckling formula, (B) our model eguation and (C) Davies's method. 





1 


2 


3 


(A) I71I 


0.0473741 


0.0969155 


0.189496 


(B) I72I 


0.047009 


0.0842059 




(C) I73I 


0.00850293 


0.0273159 


0.0492325 



In Figure [7?8l (a) we give the plot of the first critical stress values |7i| and I72I as a 
varies. In Figure [7?8l (b) and (c) we give the plot of the difference of the first critical 
stress value I72I — |7i| and the relative error ^^^|^|^^^ as a varies respectively. 

From Figure 17.81 we find that when 0.03 < a < 0.0542, the first critical stress 
value obtained by our model equations is larger than that obtained by (|7.30p . When 
a > 0.0542, the first critical stress value obtained by our model equations is smaller 
than that obtained by (|7.30p . So, according to our results, the Euler buckling formula 
gives an over-estimate of the critical buckling force when the aspet ratio is larger 
than 0.1084 and under-estimate when it is smaller than 0.1084 (although the error is 
very small, cf. Fig. 17. 8p . We also observe that when a < 0.075 the relative error is 
smaller than 5%. Also when 0.03 < a < 0.075, not only can our model equations yield 
almost the same first critical stress values with those obtained from (|7.30p . but also 
can yield corresponding buckled shapes. From the above analysis, we can say that 
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when 0.03 < a < 0.075, the Bernoulh-Euler constitutive equation induced 4th-order 
hnear ODE can match our model equations well. 

However, there are some differences between our model equations and the 4th-order 
linear ODE. Here we study a typical compressible Murnaghan material. Shearing 
is taken into account in our derivations. Moreover, our model equations only can 
yield finite number of bifurcation points when 0.03 < a < 0.0955. While, the 4th- 
order linear ODE can yield infinite countable number of bifurcation points for slender 
rectangles. 

From the above analysis, we can state that our model equation and the BernouUi- 
Euler constitutive equation induced 4th-order linear ODE are good model equations 
for the uniaxial compression of a rectangle and can match well when 0.03 < a < 
0.075. But, the 4th-order linear ODE is not a good model equation for the uniaxial 
compression of a rectangle when a > 0.075, i.e., we have determined an estimated 
interval 0.06 < 2a < 0.15 for the aspect ratio of the rectangle that the Euler's formula 
is valid. Then we believe that the Euler's formula is valid when the aspect ratio 2a is 
less than 0.15. In [25', Levinson investigated the plane strain problem of a rectangular 
parallelepiped composed of incompressible neo-Hookean material. It is found that his 
numerical estimate for buckling load can match the Euler's buckling formula when 
the aspect ratio is smaller than 0.167. 

8. Conclusions. For chosen material constants, we have found two types of 
instability phenomena in the compressions of a clamped rectangle: barrelling and 
buckling. The barrelling instability leads to a corner-like profile on the lateral bound- 
aries of the rectangle. This occurs only when the aspect ratio is relatively large. 
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i.e., when the rectangle is thick enough. More specificaUy, for the material constants 
chosen here the lower bound of the aspect ratio for barrelling is 0.376 (cf., the last 
paragraph of subsection 17. ip . Our results also reveal that this corner- like profile is 
caused by the coupling effect of the material nonlinearity and geometrical size of the 
rectangle. The analytic solution for the post-bifurcation corner-like profile is also 
obtained. When the rectangle is thin enough, buckling instability occurs instead of 
barrelling. More specifically, for the material constants chosen here the upper bound 
of the aspect ratio for buckling is O.f 9f (cf. the second paragraph below (|7.29p ). In 
the buckling case, we can recover the BernouUi-Euler constitutive equation induced 
classical 4th-order linear ODE. Numerical computations show that our model equa- 
tions and this classical 4th-order linear ODE can yield very close first critical stress 
value and the same buckling eigen-shapes when the aspect ratio is small. We point 
out that these two models are both good model equations when the aspect ratio is 
small. But when the aspect ratio is large engough, the 4th-order linear ODE cannot 
model the compression of a rectangle. 

The discovery of a lower bound of the aspect ratio for barreling and a different 
upper bound for buckling under clamped ends implies that there indeed exists a 
transition region, which agrees with the experimental results of Beatty and Hook [5] 
and Beatty and Dadras [10]. The present analytical study appears to shed certain 
light on this more than 40-years old mist. 

Appendix. 

ai ^l + r]i+ r]4, a2 = ?/o + ?/i + ?/4, as = 1 - 770 + 3ryi + 2772 + m, 

17 5 

a4 = 1 - 2r/o + 2771 + 2r]2, as = 3 -f 6771 + 2772 + 4774, ae = - + -771 + 772 + -774, 
ar = 2771 + 772 + 774, as = 7/0 + 2771 + 2774, ag = i - 770 + 771 + 774, 



aio = Vi+ V2, an = 



2 

1 + '/o + Vi + Vi - 2770 (1 + 2771 + 772 + 7^4) 



ai2 = 



ai3 



ai4 



ai5 



ai6 



ai7 



ai8 = 



ai9 



1 + 77^ + 771 + 774 - 2770 (1 + 2771 + 772 + 774) 

i + Vi+V4-Vo{i + 6771 + 2772 + 4774) 

l + 2r]^ + rii- 770 (1 + 2771 + 2772) + ?74 

277^ 2 (1 + ?7i + 774) + 77g (3 + IO771 -I- 2772 -I- 8774) - 770 (5 + IO771 + 27/2 + 8774) 

6772 

277^ - 4 (1 7^1 + 7^4) - r]^ (3 + 8771 + 4772 + 4774) + yo (7 + 14?;i + 4772 + IO774) 

6772 

(-1 + 770) (-1 - 771 - 774 + 770 (1 + 3771 + 7/2 + 2774)) 

1 - 2773 + 771 - 770 (1 + 2771 + 2772) + 774 - 4772 (771 + 774) 
677§ 

1 + 277g + 771 + 774 + 772 (2 + 4771 + 4774) - ?7o (3 H- 6771 + 2772 + 4774) 

6772 
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1 - 2?7g + ?7i + rtl (1 - 4?7i - Arfj) - 2r]o {r]2 - m) + 

1 + 27?g + 7?i + 7?4 - 27?o (2 + 4771 + 772 + 37?4) + ??o (1 + 47?1 + 4?74) 

9 + 2r7g + IO771 + 10?74 + ??g (7 + I6771 + 4772 + I2774) - 2?7o (9 + I9771 + 6772 + 18774) 

24770^ 

9 + IO771 + IO774 + ??o (9 + 12771 + 4?72 + 8774) - 2770 (7 + 15771 + 6772 + 9774) 

24772 

9 - 477g + IO771 + 7/g (7 + 4771 + 47/2) + lOr/4 - 2770 (6 + ISt/i + 6772 + 7m) 

9 + 277^ + IO771 + IO774 + ??g (5 + 16?7i + 4772 + 12774) - 27?o (10 + 21?7i + 6772 + I5774) 

247,2 

-1 + 77g - ?7i + 770 (1 + 2771 + 2772) - T?4 + 77g (-1 + 2771 + 2774) 

1 + ??o + ??i + ??4 + ?7o (1 + 27?i + 2774) - ?7o (3 + 6771 + 2772 + 4774) 

12772 

^-ril+Vi- 2% {m - Vi) + ??4 - 277g (-1 + 771 + 774) 
12772 

1 + + ^71 + ^74 + 2772 (771 + 7/4) - 2770 (2 + 4?7i + 772 + 3774) 

12r7o2 

3771 „ 3774 
1 - % + - 3770771 - 770772 + — - 2r?o??4, 

, 2 vi m 

1 3770 Vl „ ^4 „ 

— 2 + ^ - y + 3%??i + V0V2 - y + 2??o??4, 
- 2 + y ~ % - Y + '^o??! + V0V2 - -J, 

13 11 8770 277^ 27,1 IG//0A/1 4779772 2?74 

— 5- + 9771 h 3772 h 6?74 4770774, 

3 07,0 3 3 770 3 3 7,0 

11 7770 , 25771 27,1 , „ 4770772 , 167,4 2774 8770774 

4-7 ^ + 4770771 + 3772 h — —, 

6770 3 3 770 3 3 7,0 3 

8 11 7770 , 477^ 17771 27,1 4770771 4770772 , 8774 2774 

5- + ^ + 5— + 3?72 ^ + ^ , 

3 07,0 3 3 3 770 3 3 3 7,0 

16 11 5770 27,g 27,1 I6770771 477o7?2^„ 27,4 

— — + II771 — + 3772 — +8774 4770774, 

3 6770 3 3 770 3 3 770 

1 1 ^ ^ Hi ^ 47/07/1 _^ _ ^ _^ 4770774 

3 6770 3 3 3 3 33' 

4 1 277^ 77,1 4770771 4774 4770774 
3-6;^-^°- — + — -^+'^^ + ^-^' 

1 , 277§ 771 4770771 4774 , 4770774 

1 770 277^ 117,1 4770771 8774 4770774 

^"6;^"y"^ + ^ 3-+^^ + ^"^' 
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Q;42 = 1 - TT V0 + V1+ V2, a43 = 2 - + 3r]i + r/2 + 2r]4, 

2?7o 2r/o 



Q!44 



77o(-3 + 7??o-4r?^) 



3 - 8770 + 6r?^ 

Q!45 = -1 - 2?7i - 2j74 - 277^(1 + 6771 + 2r]2 + 4r?4) + ?7o(3 + 8771 + 8774) 
1 



a46 



■ (?7o(4877^ - 27(1 + 2f]i + 2774) + 16%^(-12 + 9r/i + ??2 + 8774) + 277?o(5 + 12?7i + 12774) 



(3 - 8r;o + 6r;2)3 

-2477^(9 + 32771 + 32?74) - 477^(-52 + 13577i + 3772 + I32774) + r?g(44 + 9OO771 + 9OO774)) , 

= to o^^a 2^2 (w(-36(l + 27?i + 2r?4) + 487/^(1 + 6771 + 2772 + 4r?4) 

(3 - 8770 + 6?75)^ V 

-877^(31 + 129??i + 35772 + 947?4) + 377o(67 + 154?7i + 127/2 + I42774) + 
277^(239 + 778771 + 166772 + 612774) - 77^(443 + II96771 + 1807?2 + 1016774))) , 
a48 = -1 + 2?7i + 2774 + ??o (-1 - 8771 - 8r]4) + 7?g (2 + 12r/i + 4r/2 + 87/4) . 



Di = _ (3 + 6m + 6r]4 + 47?g(377i + 772 + 2774) - 3r/o(l + 4771 + 47/4) 



^2 = o / ]^ . (-8vt + ^70(15 + 27?i + 6772 - 4774) - 27?^(3 + 7?i + 47/2 - 3774) 

377o(-l + rioj\ 

+J?o(l +V1+ V4) - 2(1 + 27/1 + 27/4))) , 



^1 = 24(_i + ^q)^3 i-H + 2 (1 + ??i + 774) + 2??2 (4 + 7??i + 3??2 + 4774) 
-??o(7+12?7i+ 4772 + 8??4)), 

O2 = (_1 + ^q)^3 (4^0 - 2 (1 + 771 + 774) + 27?^ (5 + 107?i + 27?2 + 87/4) 

+7?o (13 + 18m + 4??2 + 14774) - 2772 (11 + 2I771 + 7r/2 + I4774) ), 
^3 = Y^^Zn:^ ( - ^'^o + '?o (3 + 2r/i + 2,72) + 2 (1 + 771 + 774) 
-4770(1 + 2771+772 + 774)), 

O4 = i2(_l + ^o)^3 (2'?0 + 2 (1 + 771 + 774) - 477^ (1 + 2771 + 7?2 + m) 

+5r]l (3 + 6??i + 2772 + 4774) - 2770 (6 + 97?i + 2772 + 7774) ) , 
(?7o (1 + 4771 + 4772) - 2 (1 + 771 + 774) + 2?7^ (3771 +772 + 2774)) 

PR = 



24 (-1 + 770)773 
1 + 774) - 477^ (1 
+2?7^ (5 + ll?7i + 3772 + 8774) - ??o (11 + I8771 + 4?72 + 14774) ) , 



= 24 (-1 + 770) 773 ^^ (1 + »7i + V4) - 4770^ (1 + 3?7i + r?2 + 2774) 



_ (277g + 4 (1 + 7?i + 774) - VP (7 + 14?7i + 6772 + 8774)) 
48(-l + 77o)J7o' 
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(-1 + 4r]l - 2r]i + 2r]2 - 47?4 + 8% {r]i + 7/4)) 
24 (-1 + %) 770 

(8?7^ + 4 (1 + r?i + r?4) + Srj^ (1 + 2r/i + 2774) - 5% (3 + 6771 + 2772 + 4774)) 



48(-l + 7?o)r?2 
^ (47/g + 4 (1 + 7?i + 774) - m (5 + 10?7i + 67/2 + 4?74)) 
48 (-1 + 770) 772 
(-5 + 477g - 10771 + 27?2 - 12774 + 7?o (-4 + 8771 + 8774)) 



711 



712 



24 (-1 + 770)770 
( - 8?7^ - 4 (1 + 771 + 774) - 277^ (1 + 8771 + 8774) 



48 (-1 + 770)7,2 
+?7o (21 + 42771 + 10?72 + 32774) ) . 



1 / 9 
Hi = -IQUxUxyVx {2vx + z/4) + WyUyyVx {2ui + v^) + -U^Vxy {2vi + Ui) 

+\u^VxY (21/1 + y^) + \v^VxY {2^1 + Vi) + ^VxyV^ {2ui + 1/4) 

+WxVyVyy {2yi + Vi) + WyyV^ (A + 2vi + v^) + WyVxVxy (A + 2i/i + !y4) 
+6f/xryxVV (m + 2z/i + 1/4) + QUxVxyVy {fi + 2i/i + 1^4) 

+6UxVxVyy (a* + 2i/i + 1^4) + 9U^Uyy (A + 2/i + 2i/i + 1^4) 
+6UxyUyVy (4i/i + 21^2 + 1^4) + 6/7x1/^^14- (4i/i + 2p2 + i/4) 

+6f/xC/rl4'r (4i/i + 2i/2 + 2/4) + QUxUxyUy (^A + 2/i + 7z/i + 2i/2 + ^ 

+3U]iUYY (^A + 2/i + 7:/i + 2j^2 + + 3[/yFV^ (^A + 2/i + Ti^i + 2v2 + ^ 

+6UyVyVyy (|a + 2/i + 7z/i + 2u2 + ^ + ^ (i8U'^Vxy {i^i + 1^2) 

+36U'xfc'^V (2/1 + 2/2) + 18VxrV'r («^i + 2^2) + BC/^Vxrl^^ (A + 4z/i + 41/2) 
+3UxxV^ (A + 4!/i + 4u2) + 9Z7x?7xrVx (2z/i + 1/4) + QfeC/rVx (2z/i + 2/4) 

+9f/xC^y Vxx (2i/i + 1/4) + QUyVxVxy (/* + 2i/i + 1^4) + et/jfyVxVr (/t + 2i/i + 1/4) 
+6[/y VxxVV (/i + 2i/i + 7/4) + SU^VxY (4i/i + 2j/2 + 1^4) + SV^Vxy (4z/i + 2u2 + 1/4) 
+6t/xFt^yl^y (4i/i + 2z/2 + 1/4) + QVxVxxVy (4j/i + 2j/2 + 1^4) 

5^/4 \ , „rr Tr2 f \ . , , « , 5z^4 



+6UxUxyUy ^A + 2/i + 7i/i + 2z/2 + -j- j + ^UxxU^ ^A + 2/i + li^i +2v2+ ^ 
+3UxxV^ (|a + 2/i + 7i/i + 2i/2 + + WxVxVxx (^A + 2/i + 7z/i + 2i/2 + ^ 
+W'^xUxx (A + 2/i + 12i/i + 4i/2 + 4i/4) ) , 



if2 = i {^UlUxY (2i/i + z/4) + QUxUxxUy {2ui + 1/4) + ^UxyU^ {2vx + 2/4) 
+^C/xr'K| {2yi + z/4) + 9C/yyxVxx (2z/i + 1/4) + gC/yVxyl^Y {2ui + 1/4) 
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+ i^UxyV^ {2vi + vi) + QUxyUyVx (A + 2vi + 1/4) + iU^Vxx (A + 2vi + v^) 
+6UxUyVxy (m + 2i/i + vi) + 6UxUxyVy {fi + 2vi + v^) + WxxUyVy (m + 2vi + V4) 
+W^Vxx (A + 2/x + 2vi + vi) + WxVxVxY (4t^i + 2v2 + ^4) 
+6C/xxVxVV (4z^i + 2v2 + Vi) + WxVxxVy (4i/i + 21^2 + Vi) 

+6UxUxxVx (^X + 2ii + Ivi + 2v2 + + W^Yxx [\ + 2^ + Iv^ +2v2 + ^ 
+QVxVxyVy (^A + 2/x + 7z/i + 2i/2 + + 3VxxV^ (^X + 2ii + 7ui +2v2 + ^ 

+ i (iWlUxY {vi + 1^2) + IWxyV^ {vi + z/2) + ^WxVyVyy {vi + 1^2) 

+6C/xC/jfy^y (A + Avi + 4z^2) + ZU^Vyy (A + 4i^i + AV2) + 9[/yy V3f W (2i^i + U4) 
+9UyVxyVy {2vi + Vi) + WyVxVyy (2i^i + v^) + WxyUyVx (a* + 2i^i + ^4) 
+6f/xC/rrVx + 2z/i + 2/4) + WxUyVxy {i^ + 2z^i + 1^4) + ^UxyUy (4i^i + 2u2 + Vi) 
+WxUyUyy (4z/i + 2v2 + 2/4) + 2,UxyV^ (4z/i + 2z/2 + V4) 

+QUxVxVxY (42/1 + 22/2 + 1^4) + QUyUyyVy (^A + 2m + Ivi + 2v2 + ^ 



+6yx VxrVV (^A + 2/i + 7z/i + 2^2 + J + ^U^Vyy (^A + 2/i + 7ui + 2u2+ ^ 
+WIVyy ( a + 2/i + 7j/i + 22/2 + ^ j + OV^Vrr (A + 2/x + 12i^i + 42/2 + 4i^4) J 



iJs = QuxUxxVy {m + m) + Uxxvl - % + 2»?i + 2ry2^ + Qu^UxyVx {m + m) 

iVl + ^4) + SUyUyyVx {TJI + ??4 ) + ^U^UyV^X + ^4) + 3^ x'VyVyy 

{Vl + V4) 

+'^yy'"l - % + ??i + ^4^ + ^UxyUyVy {2t]i + r?2 + Vi) + 2uxUyyVy {2r]i + m + V4) 

+2vxVxxVy {2r}i + r?2 + m) + 2uxUyVyy {2r}i + r?2 + m) + '^u^Vxy (Sr/i + 2?72 + 774) 
3 

+ 2'^^J''^y (2^1 + 2??2 + r?4) + UxVxyVy (1 - r?o + 6r?i + 47/2 + 2r?4) 

2 / 77/1 5n4 \ 2 '''ni 5774 A o 3 „ 



+UyVxVxy (1 - 770 + 4771 + 4774) + 27i (1 + 7?7i + 2?72 + 5774) 

-\uxxUl (1 + 7771 + 2772 + 57/4) + ^ulUyy (1 + 7771 + 2772 + 5774) 

1 



+ ^UxxVl (1 + 7771 + 2772 + 5774) + (1 + 7771 + 2772 + 5774) 

+ ^"yy"y (1 + '7^1 + 2??2 + 5774) + UyVyVyy (l + Tfji + 2772 + 5774) 
+uluxx + I8771 + 67/2 + 12774^ + 2uxyVxVy (7/0 + 2 (7/1 + 774)) 

+UyVxxVy (770 + 2 (771 + 774)) + UxVxVyy {tJo + 2 (771 + 774)) , 
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QUxVyVyy {tji + r]2) + u^Vyy ~ '^0 + ^Tji + 21^2 j + iUxUxxUy {rji + 774) 

+3uyVxVxx ivi + m) + ^UyyVxVy {r]i + rn) + 6uyVxyVy + rii) + SuyVxVyy {rji + 774) 

+ulvxx Q - % + m + + 2uxUyUyy {2r]i +7/2+ m) + 4u (2r?i + 7?2 + Vi) 

3 

+2uxxVxVy {2i]i + 772 + 774) + 2uxVxxVy (2771 + 7/2 + 7/4) + -uluxy (3771 + 2772 + 774) 
3 

+ 2'^xyvl (3771 + 2772 + 774) + '^x^xy^y (1 - ?7o + 67/1 + 4772 + 27/4) 

o / 7m 5?74 \ of 7771 5774 \ n / 3 „ „ \ 

+UxyUy (^+'?2 + ^l + UxyVx ("^+^2 + — 1+ V^V^^ I - + 37/1 + 3?74 1 
+UxyUyVx (1 - 7/0 + 4771 + 4774) + '^x'^xx'^x 

(1 + 7771 + 2r]2 + 5r?4) 

+ 2'"x^a:a; (1 + 7r/i + 27/2 + 57/4) + UyUyyVy (1 + 77/1 + 27/2 + 57/4) 

+2vxVxyVy (1 + 77/1 + 2772 + 5774) + ^Vxxvl (1 + 7771 + 2t]2 + 5774) 

+ (1 + 7?7i + 2r/2 + 57/4) + ^t;>j,y (1 + 7r/i + 27/2 + 67/4) 

Q + I87/1 + 6772 + 12774^ + UxUyyVx {rjo + 2 (771 + 774)) 
+2uxUyVxy {riQ + 2 (771 + 7/4)) + UxxUyVy (770 + 2 (771 + 774)) , 

i/5 = (1 + 7771 + 27/2 + 5774) 7i27;i^ + 4 (2771 + 772 + 774) U2V1U0X 

3 

+ (1 + 77/1 + 27/2 + 5774) 7X2^0^^ + 2 (3% + 2r/2 + m) Vl^Vix 

3 

+ (1 - 7/0 + 67/1 + 47/2 + 27/4) ViUoxVlx + 2 (3^1 + 2??2 + Vi) UOx'^Vix 
+ Q - % + 27/1 + 2772^ Ul^MOxx + 6 (771 + 772) ViUOxUOxx 
+ f ^ + I87/1 + 67/2 + 127/4 j MOx^WOra, 



-f?6 = 8 + 37/1 + 37/4 j Ui7i2^ 

+3 (1 + 77/1 + 27/2 + 5774) TigUi^ + 8 (1 + 7771 + 27/2 + 57/4) ■U27;i'i;2 
+6 (1 + 77/1 + 27/2 + 5774) U1V1V3 + 12 (27/1 + ??2 + V'l) U3V1U0X 
+ 16 (27/1 +112+ Vi) U2V2UOX + 12 (27/1 + ?/2 + ?/4) W1V3MOX 
+3 (1 + 77/1 + 27/2 + 57/4) Us^Ox^ + 12 (27/1 + ?/2 + ?/4) 'W2VlUla; 
+6 (1 + 77/1 + 27/2 + 5774) U2UoxUlx + 8 (27/1 + ??2 + ??4) 'Wl^^lM2x 
+4 (1 + 7?7i + 27/2 + 5774) UiUoxU2x + 12 (t/i + 774) 7i2^t^0x 

+18 (771 + 774) V1V3V0X + 6 (770 + 2 (7/1 + 7/4)) W3^^oa;^'Ox 

+4 (770 + 2 (771 + 774)) ViU2xV0x + 12 (771 + 774) UoxU2xVOx 

„ , , / 77/1 57/4 \ 

+6 (771 + r]4)UlU2Vlx +41 — + 7/2 + — 1 U1U2V1X 
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+6 {r]i + 774) ViV2Vi^ + 6 (3r?i + 2r]2 + m) ViV2Vix 

+2 (1 - r?o + Qm + 4??2 + 2r/4) V2UoxVix + 2 + 2 (771 + r/4)) V2UoxVu 

+ (1 - Tjo + 6r]i + 4r]2 + 27^4) vmixVix + 2 + 2 {rji + 7?4)) I'l'Uix't^ix 

+6 (7?i + 774) UOxMlxt^lx + 3 (37?i + 27?2 + Vi) UQxUIxVIx 

+4 Q - r?0 + r?l + Vij U2V0xVlx + 2 (1 - 7/0 + 47/1 + 4774) U2f0x^^la: 

+ (1 - 770 + 4771 + 4r/4)Mii;i^2 ^2 ( ^ + r/2 + ^ ) ^^Oa^^^ix^ 



2 ' 2 ^ 

+3 (3771 + 2772 + 774) vi^V2x + 2 (1 - 770 + 6771 + 4r/2 + 27/4) wiUox'i^2x 

+3 (3771 + 2772 + 7/4) U0a;^t'2a: + 2 (1 + Tt/i + 27/2 + 57/4) M1M2MOXX 
+4 Q - % + 2»?1 + 27/2^ V1V2UOXX + 12 (771 + 772) V2U0xU0xx 

+6 (771 + 772) V-iUixUOxx + 2 (^^ + l^^i + + 12r/4^ Wo 

+6 (771 + 774) U2VoxUoxx + 3 (7/1 + 774) UiVixUOxx 

'1 



(1 + 7771 + 2t]2 + 57/4) VoxVlxUOxx + ( ^ - % + 27/1 + 2772 1 Vi'^Ul, 



2„ 

xx 



2„. 

XX 



+6 (771 + 7/2) t^lWOxWlxx + + l^^l + + 127?4^ WOx^Ml 

+2 (770 + 2 (771 + 774)) U2V1VOXX + 6 (7/1 + 774) U2U0xV0xx 

+2 (27/1 + % + ViVixVOxx + (1 + 77/1 + 27/2 + ^Vi) UOxVlxVOxx 

+ (r/o + 2 (7/1 + 7/4)) Mivi'i^ixx + 3 (r/i + 7/4) wi«ox'?^ixx 

+2 (27/1 + r/2 + ?/4) ViVOxVlxx + (1 + 77/1 + 27/2 + ^TJi) UQxVOxVIxx, 



8 + 37/1 + 3774 j li^ + 6 (1 + 7771 + 2r/2 + 5r/4) U4wf 

+ 18 (1 + 7771 + 2772 + 57/4) U2Wi7;3 + 24 (27/1 + 7/2 + 7/4) U4WiMox 

+ 36 (27/1 + f/2 + ?74) M2l^3U0x + 6 (1 + 77/1 + 27/2 + 5?74) U4M0X 

+20 (2r/i + 7/2 + 7/4) 'U2V1U2X + 10 (1 + 77/1 + 27/2 + 87/4) U2W0xW2x 

+ 12 (7/1 + 7/4) U^Ulx + (147/1 + 47/2 + IO7/4) Ma^lx + 18 (7/1 + 7/4) Ul^Sl'la 

+ 9 (37/1 + 27/2 + '74) V1V3V1X + 3 (1 - 7/0 + 67/1 + 47/2 + 2774) 7;3W0x^^lx 

+ 6 (770 + 2 (7/1 + 7/4)) U3Uoa;"l2; + (1 " % + 67/1 + 47/2 + 27/4) t^lW2x'?^lx 

+4 (770 + 2 (7/1 + 7/4)) ViU2xVlx + 12 (771 + 774) M0xM2x'l'lx 

+3 (37/1 + 2?72 + 7/4) U0xW2x^^lx + 2(^-7/0 + 7/1+7/4) W2^^L 



+2 (1 - 7/0 + 47/1 + 4774) 7X27;L + + »?2 + ^) f ? 

9 



+ ^ (37/1 + 27/2 + 774) vlv3x + 3 (1 - 7/0 + 6771 + 47/2 + 27/4) t^lMOx^Sx 

9 



+ - (37/1 + 27/2 + 7/4) Mox'*^3x + 2 (1 + 7771 + 27/2 + 5774) Ua^Oa 
+6 - 7/0 + 27/1 + 27/2 j l^l'l'SUOxx + 18 (7/1 + 7/2) 7;37iOxMOx 
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+6 (ryi + 772) viU2xUoxx + 2 77 + 18?7i + 6772 + 12r;4 uoxU2xUoxx 



+ {l-rio + 6771 + 4772 + 27/4) viuoxuix + - (3771 + 2772 + 774) ul^uix 

+6 (r/i + 774) U2V1V0X + 2 (?7o + 2 (771 + 774)) U2UoxVox 
+6 (771 + r/4) uivivi^ + 2 (770 + 2 (?7i + 774)) UiUoxVix 

+ 2 (1 + 7771 + 2772 + 5774) ViVOxVlx + 4 (2771 + 7^2 + 774) UOxVOxVlx 

+ (770 + 2 (771 + 774)) UiViUOxx + 3 (771 + 774) UlUOxUOxx 

+2 (2771 + 772 + 774) ViVOxUOxx + (1 + 77?1 + 27?2 + 6774) UQxVOxUOxx 

+ - (1 + 7771 + 2772 + 5774) vfvoxx + 2 (2771 + 772 + 774) ViUOxVOxx 
+ ^ (1 + 7771 + 27/2 + 57/4) uI^Vqxx, 

= 4 (1 + 7771 + 2772 + 5774) ulvi + 6 + I8771 + 6t]2 + 12774^ 7;^7;3 
+8 (2771 + 772 + 774) uluox + 36 (771 + 772) V-1V3U0X 
+6 (^-Vo + 2??i + 2772^ V3ul^ + 3 (3771 + 2772 + 774) vfu2x 

+2 (1 - 770 + 6771 + 4772 + 2774) ViU0xU2x + 3 (3771 + 27?2 + 774) ul^U2x + 

18 (771 + 774) U2V1V1X + 6 (770 + 2 (771 + 774)) U2UoxVix 

+2 (1 + 7771 + 2772 + 5774) vivl^ + 4 (27?i + 7?2 + 774) Uoxvfx + 

2 (t^o + 2 (t^i + 774)) U2V1U0XX + 6 (t^i + 7^4) U2UaxUQxx 

+ 2 (2771 + 7^2 + 7^4) ViVixUOxx + (1 + 77^1 + 2772 + 5774) UOxVlxUOxx 
+ i (1 + 7771 + 2772 + 5774) vfvixx + 2 (2771 + 772 + 774) ViUOxVlxx 

+ 2 + + 2772 + 5774) M^^^i 
18 (1 + 7771 + 2772 + 5774) U2U3V1 + 12 (1 + 7771 + 2772 + 5774) U1U4V1 
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+12 (1 + 7771 + 2772 + 5774) 7t^7;2 + 18 (1 + 77/1 + 27/2 + 5'?4) U1U2V3 
+36 (I + I8771 + 6m + 12774 ) 7;i7;2W3 + 12 ( ^ + I8771 + 6772 + I2774 I vfvA 



2 

+36 (27/1 + 7/2 + 7/4) U2U3U0X: + 24 (2771 + 7/2 + 7/4) iti7i4?iox 
+ 108 (771 + 7/2) V2V3U0X + 72 (771 + 772) 7;i7;47toj; 

+12 Q - ^0 + 2771 + 2772^ V4U0X + 8 (2771 + 772 + 774) uluix 
"^^ (~2^ + »72 + uluix + 36 (771 + 772) V1V3U1X 

+9 (3771 + 2772 + 774) V1V3U1X + 12 Q - 770 + 2771 + 2772^ vsUoxUix 

+3 (1 - 770 + 6771 + 4772 + 2774) vaUoxUix + 4 (2?7i + 772 + 774) uiW2W2x 

„ / 7771 5774 \ , „ , , 

+8 I — + 772 + — I U1U2U2X + 12 (771 + 772) V1V2U2X 

+ 12 (3771 + 2772 + 7/4) 7;i7;27t2x + 4 ( i - 770 + 2771 + 2772 ) V2U0xU2x 



^2 

+4 (1 - 770 + 6771 + 47/2 + 27/4) V2U0xU2x 

+ 3 (1 - 770 + 67/1 + 4772 + 27/4) Ui7iij;7i2x + 9 (3771 + 2772 + 7/4) UoxUixU2x 

9 

+ 2 (^'^l + + %) t^iWSx + 3 (1 - 770 + 6771 + 4772 + 2774) 7;i?i0xW3x 

9 

+ 2 (^^1 + + %) Wqx^Sx + 36 (771 + 774) U4V1V0X 

+54 (771 + 774) U2V3V0X + 12 (770 + 2 (?7i + 774)) 'U4Wox^^Ox 

+4 (1 - 770 + 4771 + 4774) Ii27i2a;7'0i; + 2 (770 + 2 (?7i + 774)) U2U2xV0x 

+36 (7/1 + 774) 7i37;i7;ia; + 48 (7/1 + 7/4) U2V2V1X 

+36 (771 + 774) U1V3V1X + 12 (7/0 + 2 (771 + 774)) usUoxVix 

+2 (1 - 770 + 4771 + 47/4) U277la;7;ia; + 6 (t/q + 2 (771 + 774)) U2UixVix 
+2 (1 - 770 + 4771 + 4774) «lU2x^^lx + 2 (770 + 2 (771 + 774)) UiU2xVlx 
+ 12 (1 + 7771 + 2772 + 5774) VsVoxVlx + 4 (2771 + 772 + 774) 'U2x^^0x1^1x 

_ fir]! 5774 \ 

+4 I — + 772 + — I U2xV0xVlx 

+5 (1 + 7771 + 2772 + 5774) U2wL + 4 (2771 + 772 + 774) uixvl^ 
7r]i ^V4\ 2 

— + 7/2 + — I Mlxflx + 30 (771 + 774) U2V1V2X 
+ 10 (770 + 2 (771 + 774)) U2UqxV2x + 6 (1 + 77/1 + 27/2 + 8774) VlVlxV2x 

+ 12 (2771 + 772 + 774) uoxVixV2x + 18 (771 + 774) 7ii7;i7;3j; 

+ 6 (7/0 + 2 (771 + 774)) ■Ui77oa;7;3a; + 6 (1 + 777i + 27/2 + 5774) Vi7;oxV3x 
+ 12 (27?i + 7?2 + 774) UOx^Ox^^Sx + 3 (?7o + 2 (?7i + 774)) WsWlWOxx 

+4 (770 + 2 (?7i + 774)) U2V2U0XX + 3 (770 + 2 (771 + 774)) U1V3U0XX 

+ 9 (7/1 + 7/4) 7i3W0xW0xx + 6 (771 + 774) U2UixUqxx + 3 (7/1 + 774) Ulti2xW0x: 
+ 6 (2771 + 7/2 + 7/4) V3V0XUOXX + (1 + 7771 + 2772 + 5774) 7i2rrW0xM0xx 
+4 (27/1 + 772 + 774) V2VixUQxx + (1 + 7771 + 2772 + 57/4) Ulx-^lxWOxx 



BIFURCATIONS OF COMPRESSIONS OF A RECTANGLE 



+2 (2?7i + ri2 + rji) ViV2xUoxx + (1 + + 2r]2 + 5t]4) uoxV2xUoxx 

+2 {r]o + 2 (?7i + r?4)) U2ViUixx + 6 (tJi + 7^4) W2W0a:Wlxa: 

+2 (2771 + ri2 + Ju) viVixUixx + (1 + 7??i + 2r/2 + 5774) uqxVixUixx 
+ + 2 (?7i + rji)) uiViU2xx + 3 (?7i + 774) uiU(3xU2xx 

+ 2 (2?7i + 772 + 774) ViVoxU2xx + (1 + 7771 + 2772 + 5774) UQxV0xU2xx 

+4 Q - % + m + '?4^ M2^0xx + 3 (1 + 7771 + 2772 + 5774) ViV^VQxx 

+6 (2771 + 772 + 774) V^UQxVOxx + 2 (2771 + 772 + 774) ViU2xVQxx 
+ (1 + 7771 + 2772 + 5774) U0xU2xV0xx + 6 (771 + 774) U2VixV0xx 

+ + 3771 + 3774 j vl^voxx + 4 Q - 770 + 771 + 774^ U1U2V1XX 

+2 (1 + 7771 + 2772 + 5774) V1V2V1XX + 4 (2771 + 772 + 774) V2U0xVi 
+2 (2771 + 772 + 774) ViUixVixx + (1 + 7771 + 2772 + 5774) UOxUixVl 
+ 6 (t^i + 774) U2V0xVlxx + 3 (771 + 774) UiVixVixx 

+2 Q + 3771 + 3774^ VOxVlxVlxx + ^ (1 + 7771 + 2772 + 5774) vf 

+2 (2?7i + 772 + 774) viUoxV2xx + ^ (1 + 7?7i + 2ri2 + 6774) ul^V2xx, 

(3 3 \ 

1 - % + - 3770771 - 770772 + - 2%??4 j ViUi 

+ {-vo - m - m) uivix + (-1 - 771 - 774) t^ox^ix 

+ l"T " y " yJ + V" 2 " y " y j 

(1 37)0 771 „ ^4 „ A 

/ 1 % 2 '71 '?4\ 
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